Interacting Preformed Cooper Pairs in Resonant Fermi Gases 
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We consider the normal phase of a strongly interacting Fermi gas, which can have either an equal 
or an unequal number of atoms in its two accessible spin states. Due to the unitarity-limited at- 
tractive interaction between particles with different spin, noncondensed Cooper pairs are formed. 
The starting point in treating preformed pairs is the Nozieres-Schmitt-Rink (NSR) theory, which 
approximates the pairs as being noninteracting. Here, we consider the effects of the interactions 
between the Cooper pairs in a Wilsonian renormalization-group scheme. Starting from the exact 
bosonic action for the pairs, we calculate the Cooper-pair self-energy by combining the NSR for- 
malism with the Wilsonian approach. We compare our findings with the recent experiments by 
Harikoshi et al. [Science 327, 442 (2010)] and Nascimbene et al. [Nature 463, 1057 (2010)], and 
find very good agreement. We also make predictions for the population-imbalanced case, that can 
be tested in experiments. 

PACS numbers: 03.75.-b, 67.40.-w, 39.25.+k 



I. INTRODUCTION 

Pair formation is a fundamental process in nature. In 
particular, it is the underlying mechanism for supercon- 
ductivity and superfluidity in interacting fermionic sys- 
tems. The impressive amount of control currently achiev- 
able in experiments with ultracold atomic quantum gases 
makes these systems ideally suited for detailed studies of 
strong correlations in many-body systems. Examples of 
such control include the cooling of the gases down to the 
nanoKelvin regime, the tuning of the interatomic interac- 
tion strength by means of Feshbach resonances, the ma- 
nipulation of the number of atoms in a particular quan- 
tum state, and the shaping of the confining potential 0- 
4]. Due to their unique properties, ultracold atomic gases 
are sometimes referred to as 'ideal quantum simulators'. 

By varying the interaction strength between fermionic 
atoms in a different internal state, it is possible to per- 
form detailed experimental studies of the continuous 
crossover between a Bardeen-Cooper-Schrieffer (BCS) 
superfluid of loosely bound Cooper pairs and a Bose- 
Einstein condensate (BEC) of tightly bound molecules 
, which results in a unified view of these two appar- 
ently different limits of superfluidity. In the intermediate 
regime of the crossover, where the scattering length of 
the interaction diverges, a novel superfluid was realized. 
In this so-called unitarity limit, the size of the Cooper 
pairs is comparable to the interparticle distance and the 
superfluid has remarkable universal properties [IlT - fl5l |. 

An early theoretical exploration of the BEC-BCS 
crossover by Leggett was performed at zero temperature 
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and made use of the mean-field BCS Ansatz 16]. The 
BCS theory takes into account a condensate of pairs, but 
ignores the noncondensed pairs. Therefore, it is not even 
suitable on a qualitative level to describe the critical tem- 
perature curve of the crossover, since at the BEC side 
superfluidity is lost due to pairs being thermally excited 
into nonzero momentum states. 

The Nozieres-Schmitt-Rink (NSR) theory for the nor- 
mal state of a balanced Fermi gas [13], which takes into 
account a noninteracting gas of noncondensed Cooper 
pairs, improves the BCS theory significantly and pro- 
vides a remarkably good description of the critical tem- 
perature T c for all interaction strengths [18|. The NSR 
theory does not change the linearized gap equation for 
the superfluid transition temperature, so that the rela- 
tionship between the (average) fermionic chemical poten- 
tial fi and T c is unchanged with respect to BCS theory. 
However, the NSR theory does change the equation of 
state, which affects the total particle density n(fi, T). As 
a result, the ratio fceTc/eF is altered, with the Fermi en- 
ergy given by ep = ?i 2 (37r 2 n) 2 / 3 /2m and m the mass of 
the fermions. At unitarity, NSR theory predicts rather 
large values for the ratios k^T c / /i and ef/jU compared to 
Monte-Carlo (MC) results. Nevertheless, the NSR pre- 
diction fceTc = 0. 23eF is still quite close to the MC result 
of kftT c = 0.15ep [19(. It was also shown that away from 
the critical temperature the NSR theory gives excellent 
agreement with accurate thermodynamic measurements 
on the unitary balanced Fermi gas [ill - flol [20l l2lT ] . 

In the beginning of 2006, two experimental groups 
performed the first experiments on an ultracold atomic 
Fermi gas with a po pul ation imbalance in its two acces- 
sible spin states [22l . 23] . The phase diagram was exper- 
imentally shown to be dominated by a tricritical point 
[2~i| . that was predicted from mean- field theory, NSR 
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theory and renormalization-group theory (25l - [27j . How- 
ever, in Ref. [27| it was already noticed that the NSR 
theory breaks down for small spin imbalances. Namely, 
near the critical temperature, the NSR theory predicts 
a negative polarization p = (n+ — u_)/(n+ — nJ) for a 
positive chemical potential difference —//_), which 
corresponds to a compressibility matrix —d 2 uj/d[j, a d[i a -i 
that is not positive definite [271 ]. Here, n± are the atomic 
densities for the two spin states a = ±, \i a is the chemi- 
cal potential for spin state a and lj is the thermodynamic 
potential density of the Fermi gas. It is a quite unsatis- 
factory situation that the NSR theory, which gives such 
a good agreement with accurate thermodynamic exper- 
iments for the strongly- interacting normal state of the 
Fermi gas [H], EH l2lj |. already gives unphysical results 
for even the smallest population imbalances. 

For this reason, we improve the theory of Cooper- 
pairs by taking also into account the effect of the in- 
teractions between the noncondensed pairs. If the mi- 
croscopic fermionic action is exactly transformed into a 
Cooper-pair action by means of the so-called Hubbard- 
Stratonovich transformation [28|, [29[, then the result- 
ing action not only contains a noninteracting part, but 
also two-pair interactions, three-pair interactions, and 
all higher-order interactions. In the BEC limit of the 
crossover, the tightly bound pairs interact repulsively 
with a scattering length given by 0.6a with a the scatter- 
ing length of the interaction between the fermions (30j . 
A Popov theory for the composite bosonic pairs that 
includes the pair interaction effects was formulated by 
Pieri and Strinati, leading in the BEC regime to Popov's 
results for point-like bosons [3l], H2|- Below the criti- 
cal temperature, also the Bogoliubov theory for inter- 
acting Cooper pairs was studied [2(1 |33|, Other 
strong-coupling approaches that go beyond the NSR the- 
ory include so-called self-consistent ladder approxima- 
tions [H, HI] and Monte-Carlo calculations 
It is somewhat surprising that these calculations do not 
seem to lead to a better agreement with most recent ex- 
periments on the equation of state for a unitary Fermi 
gas than the NSR theory [H[2l[. 

Viewing the normal state of a strongly interacting 
Fermi mixture as a gas of interacting Cooper pairs 
is complementary to a Fermi- liquid picture [41] . The 
Fermi-liquid picture focusses on the single-particle corre- 
lation function, or fermionic Green's function, while the 
Cooper-liquid picture focusses on the two-particle corre- 
lation function, or the pair Green's function. Possible 
differences between the two pictures arise only from dif- 
ferent approximation schemes. The additional advantage 
of the Cooper-liquid picture is that it gives a clear con- 
dition for the transition to the superfluid phase, namely 
when the (effective or renormalized) chemical potential 
of the Cooper-pairs goes to zero. 

To study interaction effects in an interacting Bose gas 
in a nonperturbative manner, the renormalization group 
(RG) approach is an established approach [42j]. RG stud- 
ies for interacting bosons on the verge of becoming su- 



perfluid have increased our understanding of the result- 
ing phase transition [43l - |47j . In this article, we perform a 
renormalization-group study of interacting Cooper pairs 
in the unitarity limit. To this end, we generalize the 
Wilsonian RG theory for point-like bosons to the more 
complicated case of Cooper pairs. The article is orga- 
nized as follows. First, we briefly discuss the exact action 
for the Cooper pairs that can be derived from the micro- 
scopic fermionic action. Then, we set up the Wilsonian 
renormalization scheme for Cooper pairs in order to cal- 
culate the effects of the Cooper-pair interactions. In par- 
ticular, we calculate the self-energy of the Cooper pairs. 
We also compare our results with the recent experiments 
by Harikoshi et al. [3| and Nascimbene et al. [l5| . The 
agreement with these detailed experiments on the equa- 
tion of state of the unitary Fermi gas turns out to be 
very good. We also make predictions for the equation of 
state of the imbalanced gas, which can be experimentally 
tested in the near future. 



II. EXACT ACTION FOR COOPER PAIRS 

We start from the microscopic action for an ultracold 
two-component Fermi gas with a local interaction 

sy>*,<l>] = -£ /^^(O^-^OMO 
+v J d$ (!) 

where 0o-(£) is the fermionic field associated with the an- 
nihilation of a particle with spin a at £ = (t, x). Here, 
x denotes the position and r the imaginary time. More- 
over, Vo is the strength of the local interaction, while the 
Fourier transform of the inverse noninteracting Green's 
function G®^ 1 is given by ?iG .^ k = ihuj n — e k + /J<x, 
with fi a the fermionic chemical potential for spin state 
<7, e k = h 2 'k 2 /2m the kinetic energy, and uj n the odd 
fermionic Matsubara frequencies. Note that in Eq. (fTJ) 
the integration over r is from to hjk^T with T the 
temperature. For atomic gases the spin label a refers to 
two hyperfine states that are used in the experiment to 
realize the two-component gas. The action of Eq. (fT]) was 
previously used as the starting point for a RG study of 
the population-imbalanced Fermi gas in its normal phase 

For an attractive interaction, the purely fermionic mi- 
croscopic action can be exactly transformed into a Bosc- 
Fermi action that contains the pairing field A(x, r) by 
means of the Hubbard-Stratonovich transformation. For 
the precise procedure, see e.g. Ref. Q. The result is 

,0,A*,A] = (2) 

-£ / C(€)^ _1 (€.0^(0 

a=± J 

+ J ^{-^|^ + A*(0<MO<MO + c.cJ, 
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where the last two terms of the second line indeed show 
that two fermions of opposite spin can form a pair, or that 
a pair can decay into two fermions. The action of Eq. ([2]) 
can be interpreted as an interacting Bose-Fermi mixture 
[48l |. and has been the starting point of the functional 
renormalization group studies in Refs. [49l l5(il|. 

Here, we follow a different route. We start by integrat- 
ing out the fermions exactly, resulting in an exact micro- 
scopic action for the Cooper pairs. We call this action for 
the pairs microscopic, because no pair fluctuation effects 
have been taken into account yet. Incorporating these 
fluctuations exactly would consequently result in the ex- 
act effective action for the pairs. The microscopic Cooper 
pair action has a very rich structure in momentum and 
frequency space, i.e., a highly nonlocal character, which 
in particular is true for the obtained Cooper-pair vertices. 
In this study, we take the complicated structure of the 
microscopic Cooper-pair propagator fully into account 
when performing the RG calculations. However, when 
calculating the effective pair propagator with our RG, we 
only consider the renormalization of the momentum- and 
frequency- independent part of the pair propagator, i.e., 
we do not consider the renormalization of the pair effec- 
tive mass. After the exact integration over the fermionic 
fields, Eq. ^ results in (see again e.g. Ref. 1) 

S[A*,A] = - f dfJ^P^ -OTrlogt-G- 1 ), (3) 
J Vq 

where the Nambu space inverse Green's function G _1 
is given by G n X = G+~ 1 (£,£')i G£ = -GPt 1 
hG^ 1 = -A(£)<S(£ - £') and G^ 1 = G^*. By expand- 
ing the logarithm in Eq. ([3J in powers of A, we obtain an 
infinite series that prohibits an exact solution to the prob- 
lem, so that approximations have to be made in order to 
proceed. By performing the mean-field, or saddle-point, 
approximation, the full path integral over the bosonic 
field A(£) is simply approximated by the value of the in- 
tegrand associated with the global minimum S*[Aq, Aq]. 
This approximation results in the well-known BCS ther- 
modynamic potential. 

Going beyond mean-field theory, the next step is to 
perform a Gaussian, or random-phase approximation. In 
the normal phase, this means that Eq. ([2]) is expanded 
up to second order in the pairing field, and that the re- 
sulting Gaussian functional integral is performed exactly. 
By neglecting the higher-order contributions, the Cooper 
pairs are thus physically approximated as forming a non- 
interacting gas. The resulting theory is also called the 
Nozieres-Schmitt-Rink approximation [l7j . and it has 
been applied with success to the study of thermody- 
namic properties above and below the critical temper- 
ature pol . [33l [HI • It is namely possible to generalize this 
theory also to the superfluid state by making in Eq. @ 
the substitution A(£) = A + A'(£), expanding Eq. © 
up to second order in the fluctuations A' , and performing 
the resulting Gaussian functional integral. 

In this article, we only study the normal phase. The 
inverse propagator for the noncondensed Cooper pairs 
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FIG. 1: Diagrammatic representation of a) the bare Cooper- 
pair propagator and b) the bare Cooper-pair interaction. The 
Cooper pairs are represented by thick lines, while the thin 
lines correspond to fermionic propagators. 



G^ that follows from the quadratic part in the pairing 
field of Eq. ([3]), is given by (see e.g. Ref. @) 

%G\- X {iuj n ,k)= (4) 

m , 1 W l-/+(ek')-/-(ek-k») 1 ] 
AnTi 2 a V ^-f \ -%huj n + e k ' + £k-k' - 2e k ' J ' 

with V the volume, a the scattering length, fi = 
(/j+ + /U_)/2 the avarage chemical potential, f a (e) — 
l/je^ 6- ' 1 ^ + 1} the Fermi distribution function, j3 — 
l/k^T, and uj„ an even bosonic Matsubara frequency. 
The Feynman diagram that corresponds to the Cooper- 
pair propagator is shown in Fig. [TJi. We call G^ _1 
the bare or microscopic propagator, indicating that no 
Cooper-pair interaction effects have been taken into ac- 
count yet. Note that the bare propagator is exact, in the 
sense that it follows from an exact transformation of the 
fermionic action. With our RG approach we can con- 
sequently systematically include Cooper-pair interaction 
effects that lead to self-energy corrections to the bare 
propagator. 

The Cooper-pair interaction V® follows from the quar- 
tic part in the pairing field of Eq. ([3j and is diagrammat- 
ically represented in Fig. [Tb [3l|, |32[. Here, we do not 
take the full frequency and momentum dependence of 
the Cooper-pair interaction vertex into account, but we 
consider only two external frequencies and momenta to 
be nonzero, namely either u)\ — —u>2 and ki = — k2, or 
CJ3 = — CJ4 and k3 = — k 4 , where the labeling is given 
in Fig. [T]d. This specific choice corresponds physically 
to considering only zero center-of-mass frequencies and 
momenta, which is motivated in the next section. The 
resulting expression is given by 

Vi(iuj n , k) = V£G v (iu n , k) = x (5) 

Ur -,n',k' t - T + : -n',-k' LT -,n'+n,k'+k l - T + ,-ri'-n,-k'-ki 

n' ,k' 

where we have defined = V^(0, 0), so that Gy en- 
capsulates the considered (relative) momentum and fre- 
quency dependence of the Cooper-pair interaction. The 



4 



Matsubara sum over odd fermionic frequencies n' in 
Eq. (O is readily performed analytically, but results in 
a somewhat cumbersome expression. We call V®(iuj n , k) 
the bare or microscopic interaction, in order to make the 
distinction with the effective or renormalized Cooper-pair 
interaction, which includes the effect of Cooper-pair fluc- 
tuations and that is calculated in the next section during 
the RG flow. 

Due to the repulsive interaction between the Cooper 
pairs, they acquire a self-energy £a, of which the 
momentum- and frequency-independent part is most rel- 
evant for the RG flow. So the full propagator becomes 
G^ 1 (iw„,fc) = G A _1 (ia; n ,fc) — £a- We define the bare 
Cooper-pair chemical potential as /j, a = 7iG A -1 (0,0), 
while the renormalized chemical potential is given by 
A'A = /^a — ^Ea- The renormalized chemical potential 
thus includes the self-energy effects. As a result, the 
full Cooper-pair propagator Ga depends on the renor- 
malized chemical potential /ia- In this study, we do not 
take the frequency and momentum dependence of the 
Cooper-pair self-energy into account. Doing this would 
result also in a renormalization of the effective mass of 
the Cooper pairs, an interesting topic for further study. 
Note that the Cooper-pairing fields do not have the same 
units as the fields for point-like bosons, resulting also in a 
different unit for the propagator. Therefore, the present 
definitions of the Cooper-pair chemical potential and the 
self-energy do not have the unit of energy. However, we 
still think our nomenclature is appropriate, due to the 
physical and mathematical analogy with the correspond- 
ing concepts for point-like bosons. 

A complementary physical meaning for the chemical 
potential of the Cooper pairs is obtained by realizing 
that the noncondensed Cooper pairs mediate an inter- 
action between the fermions as follows from the action in 
Eq. ([2]). The bare Cooper pair propagator of Eq. ((4]) 
is indeed equivalent to the many-body transition ma- 
trix for the fermions in the ladder-diagram approxima- 
tion, which we call the bare many-body transition ma- 
trix. At zero energy and momentum, we have for the 
bare many-body scattering length a^ B that 7iG A (0, 0) = 
l//x A = T^ B (0, 0) = 4irh 2 a^ B /m. The renormalized 
chemical potential for the Cooper pairs then corresponds 
to a renormalized transition matrix for the fermions that 
includes more Feynman diagrams than the bare one, 
for which the renormalized many-body scattering length 
omb is given by fiA = m/Airh 2 aMB- 

III. RENORMALIZATION FORMALISM 

To treat the interaction effects of the Cooper pairs in 
a nonperturbative manner we use the Wilsonian renor- 
malization group approach [42| . The procedure goes as 
follows. First an integration is performed over degrees of 
freedom in a high momentum shell of infinitesimal width. 
The result of this integration is consequently absorbed 
into various coupling constants of the theory, which are 



said to flow if degrees of freedom in subsequent momen- 
tum shells are integrated out. We could also perform an- 
other step, namely a rescaling of the momenta, frequen- 
cies, and fields. This is convenient, if we would wish to 
treat universal properties of critical phenomena such as 
critical exponents by looking at so-called RG fixed points 
42]. In this article, however, we do not wish to calcu- 
late universal critical exponents, but rather quantities 
like the self-energy of the Cooper pairs, for which rescal- 
ing is not particularly useful. The renormalization group 
then serves as a nonperturbative method to iteratively 
solve a many-body problem, rather than as a mapping 
between actions with the same high-momentum cutoff 
from which critical scaling relations can be derived. 

Thus, the first step of the method is to evaluate the 
Feynman diagrams that renormalize the coupling con- 
stants of interest, while keeping the integration over 
the internal momenta restricted to the considered high- 
momentum shell. Only one-loop diagrams contribute to 
the flow, because the width of the momentum shell is in- 
finitcsimally small and each loop introduces an additional 
factor proportional to the infinitesimal width. Although 
the one-loop structure of the infinitesimal Wilsonian RG 
is exact, it does not mean that it is easy to also obtain ex- 
act results, since this would require the consideration of 
an infinite number of coupling constants. Although the 
latter is usually not possible in practice, the RG distin- 
guishes between the relevance of the coupling constants, 
so that a small set of them may already lead to accurate 
results. 

The simplest RG calculation that gives nontrivial re- 
sults treats the renormalization of the chemical potential 
\x a and the interaction strength Va . It ignores the three- 
pair interactions and higher. The flow equations for these 
two coupling constants can be derived along exacly the 
same lines as for point-like bosons, for which detailed ac- 
counts can be found in Refs. @, The corresponding 
one-loop Feynman diagrams are diagrammatically repre- 
sented in Fig. [5] The main difference with point parti- 
cles is that the frequency dependence of the Cooper-pair 
propagator is more complicated. As a result, the Mat- 
subara sums of the one-loop Feynman diagrams cannot 
be performed analytically anymore, but have to be eval- 
uated numerically within each momentum shell. 

The RG flow of the Cooper-pair interaction strength 
and chemical potential are determined by the following 
set of coupled differential equations 

= Pufan a, Va), = Pv(l,VA, V A ) , (6) 

where the '/3-functions' are given by 

Pfi = -k[^j^2^GA{i^n,k u ^A), (7) 

n 

Pv = fc i '^{S(A; i ,/iA)+4n(fc i ,/ZA)}. (8) 

Here, n and S are the so-called 'bubble' and 'ladder' con- 
tributions to the effective Cooper-pair interaction, which 
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FIG. 2: Diagrammatic representation of the l j3 functions', a) 
Feynman diagram determing the self-energy of the Cooper 
pairs, b) Feynman diagrams renormalizing the Cooper-pair 
interaction. The middle diagram is also called the ladder 
diagram, the right diagram is called the bubble diagram. Note 
that the lines are thick and correspond to the Cooper-pair 
propagator. 



are explained in the next paragraph. Moreover, k\ de- 
notes the wavevector of the Cooper pairs in the shell of 
infinitesimal width. This wavevector is parametrized by 
the flow parameter I, and we start the RG flow at the 
high- momentum cutoff TiA and decrease as ki = Ae~ l . In 
addition, k[ is the derivative of k with respect to I. Solv- 
ing Eq. ([6|) for increasing I means that we are includ- 
ing the effect of pair fluctuations with lower and lower 
momenta, while due to the coupling of the differential 
equations we automatically generate an infinite number 
of Feynman diagrams, showing the nonperturbative na- 
ture of the RG. The initial conditions for Eq. ([6]) are 
H A (l = 0) = /4 = hG A _1 (0,0) and V A (l = 0) = v£, 
which are calculated from Eqs. jl]) and (0. 

The one-loop expression for the renormalization of the 
chemical potential in Eq. (JT)) that determines the self- 
energy of the Cooper pairs, has a clear physical meaning, 
since it is seen to be proportional to the renormalized 
pair interaction strength and to the density of Cooper 
pairs. The 'bubble' diagram II(fc, /ia) describes the effect 
of particle-hole fluctuations on the effective Cooper-pair 
interaction, where these particles are now Cooper pairs. 
It is given by 



II(fc, (j, A ) = — ^ G A (iuj n , k, ha) 



(9) 



The 'ladder diagram' describes the Bose-enhanced scat- 
tering of the bosonic Cooper pairs, given by 

3(fc,jU A ) = -j ^G v (iuj n ,k) 2 \GA(ioJn,k,fiA)\ 2 , (10) 

n 

where the momentum and frequency dependence of the 
interaction in Eq. ([5]) . i.e. Gy, is seen to enter. This fre- 
quency and momentum dependence is important, since 
otherwise Eq. (fl"U| ultimately would lead to an ultravi- 
olet divergence. This divergence physically arises from 
approximating the pair interaction as a point interac- 
tion, which is therefore insufficient. We also note that 
the self-energy diagram of Fig. [5^ and Eq. ([7]) , and the 



bubble diagram of Fig. [2Jd and Eq. ([9]) do not lead to 
divergencies. As a result, our present scheme for includ- 
ing the Cooper-pair interactions is the minimal choice 
for obtaining divergence-free, or equivalently, cutoff in- 
dependent results. 

The structure of Eqs. © to © is analogous to that 
for the RG equations in the normal state of Ref. [4J], 
where point-like bosons are treated, but there are also a 
few differences. In that reference, a trivial scaling was 
performed on the chemical potential and the interac- 
tion strength, which amounts only to a direct rewrit- 
ing of the differential equations. Second: there, also 
three-body interaction effects were considered, which are 
ignored here. Third: the present Cooper-pair propa- 
gator is more complicated than the atom propagator, 
so that the Matsubara sums in Eq. ([7]) and ((U cannot 
be done analytically. This difference can be avoided 
by appoximating the inverse Cooper-pair propagator 
as G A 1 (iw n , k, (ia) = (ihuj n - h 2 \t 2 /2toa + MaV^a, 
meaning that we would perform a low-energy and long- 
wavelength expansion with toa the effective mass of the 
Cooper pairs and {i'a = ^AMA having the unit of energy. 
In this article, however, we take the full frequency and 
momentum dependence of the propagator into account, 
since this gives the most accurate results. 

The last difference with Ref. [3] involves the ladder 
diagram of Eq. (|10p . Such ladder diagrams are known to 
lead to an ultraviolet divergence, if the interaction is ap- 
proximated as a point interaction, i.e., when Gy — 1. For 
interacting bosonic or fermionic atoms, the point inter- 
action approximation is often used. The problem is then 
most easily solved by choosing a value for the microscopic 
interaction that cancels the divergence and that leads in 
the two-body limit of the theory to the scattering length 
that is obtained from experiments. For our present RG 
flow of the interacting Cooper pairs, this procedure is not 
possible, since the bare Cooper-pair interaction strength 
is calculated exactly from the transformed microscopic 
action, and is not a free parameter that we can try to 
match to some experimentally known quantity. As a re- 
sult, we really need to go beyond the point interaction 
approximation and take into account the momentum and 
frequency dependence of the pair interaction in order to 
obtain physical results |3l| . 



IV. RESULTS 

We compare our RG equations with recent detailed 
thermodynamic experiments that have measured the 
equation of state for a homogeneous unitary Fermi gas 
[Til HH . A relevant question is how our theory changes 
the results that follow from the NSR theory. As also 
found in previous studies, the results from the NSR ap- 
proximation are in excellent agreement with experiments 
on the thermodynamics of the balanced Fermi gas at 
unitarity [2l| . However, a fundamental problem with 
the NSR theory is that it gives already for small pop- 
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FIG. 3: (Color online) Bare many-body scattering length 
a MB (dashed line) and fully renormalized many-body scat- 
tering length a^e (full line) as a function of temperature 
T in the unitarity limit for the balanced Fermi gas. The 
Fermi wavevector /cf and Fermi temperature Tp are directly 
related to the total particle density through the Fermi en- 
ergy e F = k B T F = h 2 kl/2m = ft 2 (37r 2 n) 2/3 /2m, which is 
calculated with the RG approach. The many-body scattering 
length is inversely proportional to the Cooper-pair chemical 
potential, and the difference between the bare and the renor- 
malized value is due to the repulsive Cooper-pair interaction 
effects. 



ulation imbalances rise to unphysical compressibilities 
—d 2 oj/dii a dn <7 > [I?]] • Moreover, the ratio of the critical 
temperature and the chemical potential is not so accu- 
rate for the NSR approximation, since it is the same as 
for mean-field theory, giving rise to fcB^MF = 0.66// for 
the balanced Fermi gas. This follows from the observa- 
tion that when G A -1 (0,0) = Ma = 171 1 ^mb = °> 
then both mean-field theory and NSR theory predict a 
transition to the supcrfluid state, i.e., the condensation of 
Cooper pairs. The result deviates more than a factor of 
two from the Monte-Carlo result k B T^ c = 0.31/x [HI. So 
a theory that goes beyond NSR ideally should not have 
a large effect on the equation of state for the balanced 
case, should significantly lower the critical temperature, 
and give rise to physical results in the imbalanced case. 

Studying the thermodynamics of the unitary Fermi gas 
in the grand-canonical ensemble requires the knowledge 
of the thermodynamic potential. The NSR theory of the 
strongly interacting normal state gives rise to two contri- 
butions to the thermodynamic potential, namely a con- 
tribution describing an ideal gas of fermions and a contri- 
bution describing an ideal gas of noncondensed Cooper 
pairs [17} . The ideal Fermi gas contribution follows di- 
rectly from Eq. ([3]) as the part that is independent of the 
Cooper-pair field A, and is therefore also not renormal- 
ized. It is given by 



w ig 



•igCT,AO = - 



,-0(e k -^ (11) 



The contribution to the thermodynamic potential den- 
sity due to the Cooper pairs is given by the one-loop 
expression 



~~dT 



k 2 x , 



(12) 



where the first minus sign on the right-hand side is only 
present when k\ is a decreasing function. Note that 
this last expression gives precisely the differential form 
of the NSR contribution to the thermodynamic poten- 
tial density when the Cooper-pair chemical potential is 
not renormalized (/ia = Ma)j i- e -j when we consider the 
Cooper pairs to be noninteracting. To evaluate Eqs. ([6]) 
and (1121) numerically, it is convenient to perform contour 
integration, leading to the results in Eqs. (|A2|) . (|A3[) and 
(|A4p of the appendix [13] ■ If the exact Cooper-pair prop- 
agator is inserted in Eq. (|T2"j) . then the exact thermody- 
namic potential density is obtained. However, this would 
require the treatment of all n-body interactions, which is 
presently out of reach. 

In this article, we have studied the effect of the two- 
pair interactions on the thermodynamic potential for var- 
ious temperatures. This was done by simultaneously 
solving Eqs. ((6]) and (fT2")l numerically. We obtain the 
renormalized Cooper-pair chemical potential at the end 
of the flow, /iA.cxj = (J-a(1 — ► 00), which incorporates 
the pair self-energy effects. The results are shown in 
Fig. [3] for the balanced case (fi„ = fi), where the (renor- 
malized) many-body scattering length is plotted, which 
was introduced in Section |TTJ This scattering length is 
inversely proportional to the (renormalized) pair chem- 
ical potential, namely 1/fcFiMB = ^kTi 2 iik/nik-p with 
kp = (2meF) 1//2 /?i the Fermi wavevector. In Fig. [3J 
the initial or bare many-body scattering length a^ B 
and the final or fully renormalized many-body scatter- 
ing length a^ B are plotted as a function of tempera- 
ture. We see that a^ B 1 is more negative than a^^ 1 , 
which means that the Cooper-pair self-energy, given by 
?i5]a,oo = A*a — Ma, 007 is positive, caused by the repulsive 
pair interaction. Alternatively, interpreting the Cooper- 
pair chemical potential as being proportional to the in- 
verse many-body transition matrix, we have for the con- 
sidered temperature range that |a^J B | < |a^ B |, which 
means that the fluctuation and interaction effects make 
the gas less strongly interacting. 

To see whether the selfenergy of the Cooper pairs, or 
the reduced many-body scattering length, has observable 
experimental consequences, we calculate the equation of 
state, which has been accurately measured by Harikoshi 
et al. [IH and Nascimbene et al. [l5| for the balanced 
case. We have that the total thermodynamic potential 
density is given by 



u)(T, n a ) 



n(r, /i 
v 



■ +wa, 00 (T, / 0, (13) 



with Q the thermodynamic potential and uja,qc(T, fJ- a ) — 
ua(1 — > 00). From u>, all other thermodynamic quanti- 



l " 1 " r 



0.8 - 




0.4 U 1 1 1 L 

0.3 0.4 0.5 

T/T F 



FIG. 4: (Color online) Equation of state for the normal phase 
of a strongly interacting balanced Fermi gas at unitarity in the 
canonical ensemble. The energy E of the gas is calculated as 
a function of temperature T with the renormalization group 
approach (full line) and the Nozieres-Schmitt-Rink approach 
(dashed line). The squares are the experimental results of 
Harikoshi et al. [T3 |. the triangles of Nascimbene et al. [l5l |. 
and the dots are the Monte-Carlo results of Burovski et al. 
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FIG. 5: (Color online) Equation of state for the normal phase 
of a strongly interacting balanced Fermi gas at unitarity in 
the grand-canonical ensemble. The pressure P — —Q/V of 
the gas is calculated as a function of the inverse fugacity e - ^ 
with the renormalization group approach (full line) and the 
Nozieres-Schmitt-Rink approach (dashed line). The pressure 
of the ideal two-component Fermi gas is given by Pi g . The 
triangles are the experimental results of Nascimbene et al. 




ties of interest can be obtained by the standard thermo- 
dynamic relations. Of particular interest are the particle 
densities n a = —dcu/d^a- and the energy 

E = Q + fi + n + V + fi-n-V + TS, (14) 

where S is the entropy, given by S = —dfl/dT. The 
results of our calculations for the balanced case are 
shown in Fig. 2J where also the experimental results of 
Refs. 114 [loll are given, as well as the Monte-Carlo results 
of Ref~[ll7and the results from NSR theory 0, [1J|. It 
is seen that the inclusion of the Cooper-pair self-energy 
with the RG theory results in a deviation from the NSR 
theory that is comparable with the size of the error bar 
in the experiments and Monte-Carlo calculations. Fig. 0] 
plots the thermodynamic quantities for the canonical en- 
semble, while Fig.[5]shows the comparison with the data 
from Nascimbene et al. in the grand-canonical ensemble, 
which gives the most direct comparison with their mea- 
surements [Hj]. Fig.[S]shows that in the normal state, the 
agreement with the NSR theory is perfect over a very 
large temperature range. This is remarkable, since the 
NSR theory is not exact, and in the strongly- interacting 
regime deviations might have been expected. With our 
RG theory we calculate the leading pair self-energy ef- 
fects beyond the NSR calculation, and find that the ef- 
fect on the equation of state is quite small, as seen in 
Fig [5] However, the effect is still observable and the 
agreement with experiments becomes worse. A theoreti- 
cal explanation for this effect could be that, although we 
go far beyond the NSR approximation, there are more ef- 
fects that play a small quantitative role in the comparison 



with experiments, like for example the renormalization of 
the effective Cooper-pair mass and three-pair interaction 
effects. These are interesting topics for further study. 

We have also calculated the equation of state for the 
imbalanced Fermi gas. In Fig. [SJ we show the pressure 
P as a function of h/fi for the temperature T/fi — 0.75, 
where h = — //_)/2. At this temperature, the NSR 
approximation predicts a negative polarization for pos- 
itive h, which is an unphysical result. Our RG theory 
that treats the interaction effects beyond NSR theory 
does not have this problem. We see that, as a result, our 
RG theory, the NSR theory and the mean-field theory 
give very different results for the pressure of the imbal- 
anced Fermi gas in Fig. [6] This pressure has recently 
been measured at zero temperature, where good agree- 
ment with Monte-Carlo calculations was obtained (4l| . 
The pressure could also be measured above the critical 
temperature, giving rise to a sensitive test for theories of 
the imbalanced Fermi gas at nonzero temperatures. 

Finally, we briefly discuss the effect of the Cooper-pair 
interactions on the critical temperature for the balanced 
Fermi gas. In the unitarity regime, the effective two- 
pair interaction is repulsive in the normal state. Taking 
into account the repulsive two-pair interaction will lower 
the ratios k^T c / fx and cf/a 4 , because physically the re- 
pulsive interactions lower the effective chemical poten- 
tial of the noncondensed Cooper pairs. As a result, the 
density of noncondensed Cooper pairs is lowered, which 
decreases the total density and the critical temperature 
for condensation. Namely, when G^ _1 (0,0) = ^ = 
m/ Airh 2 a^ B — 0, then both mean-field theory and NSR 
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FIG. 6: (Color online) Equation of state for the normal phase 
of a strongly interacting imbalanced Fermi gas at unitarity in 
the grand-canonical ensemble. The pressure P = —u(T,jj,,h) 
of the gas is calculated at a temperature T = 0.75/1 as a func- 
tion of the normalized chemical potential difference h//j, = 
(/!+ — /!-)/ (/Lt++/x_) with the renormalization group approach 
(full line), the Nozieres-Schmitt-Rink approach (dashed line), 
and for the ideal Fermi gas (dashed-dotted line). For each 
curve, the pressure of the imbalanced gas is normalized by the 
corresponding pressure of the balanced gas Pq = —u{T, fj,, 0). 



theory predict a transition to the superfluid state, i.e., 
the condensation of Cooper pairs. This gives rise to 
^b-^mf = 0.66//, which deviates more than a factor of 
two from the Monte-Carlo result feB^MC = 0-31// fl9j - 
Upon lowering the temperature from T^ F , /z^ becomes 
positive, but due to the repulsive interactions the renor- 
malized chemical potential is lowered, and at the end of 
the flow /xa,oo might still be negative. Therefore, the 
critical condition in the presence of Cooper-pair inter- 
actions becomes /ia.oo = 0. With our present theory, 
this results in fcB?R G = 0.43/i, significantly closer to the 
Monte-Carlo results than the mean-field and the NSR 
result. Note that this result can be further improved 
by performing a RG calculation for the superfluid state. 
Namely, close to T c we then initially have /za(0 > 0, and 
the RG flow starts out in a superfluid state. Due to the 
repulsive interactions, the effective chemical potential is 
then again lowered, and the system can flow into the nor- 
mal state [44j . This calculation thus requires a superfluid 
RG, and although progress has been made in this direc- 
tion using various approximations (such as a low-energy 
expansion of the Cooper-pair propagator) , the full calcu- 
lation is much more tedious than for the normal phase 
and is prospect for future research. At this point, we can 
also explain our choice for the range of the temperature 
axis in Fig. 2J since for lower temperatures it would be 
better to perform the partly superfluid RG. 



V. DISCUSSION AND OUTLOOK 



To summarize our results, we have calculated 
the thermodynamic properties of strongly-interacting 
Fermi gases. This was achieved by combining two 
well-established strong-coupling methods, namely the 
Nozieres-Schmitt-Rink theory for the treatment of non- 
condensed Cooper-pairs together with the Wilsonian 
renormalization group scheme for treating the pair inter- 
action effects. Our resulting theory incorporates fluctua- 
tion effects far beyond the NSR theory, and the obtained 
equation of state has been compared for the balanced 
Fermi gas with the experiments of Refs. [3, [H|. The 
difference between the equation of state from our RG the- 
ory and NSR theory has about the same size as the error 
bars in the experiment of Ref. [l4j]. As a result, taking 
the theoretical step from the noninteracting NSR theory 
to the incorporation of pair interaction effects maintains 
the good agreement with accurate thermodynamic exper- 
iments. A detailed comparison between theory and the 
experiment of Ref. [l5| shows that the results from the 
NSR theory are in perfect agreement over a large temper- 
ature range, and the inclusion of self-energy effects makes 
this agreement somewhat worse. In order to try to under- 
stand this result from the theoretical side, the accuracy 
of the calculation could be further enhanced, namely by 
calculation of next-to-leading order effects, such as the 
renormalization of the Cooper-pair mass and the three- 
pair interaction, which might also have quantitative ef- 
fects. These are projects for the future. 

The NSR approximation fails in describing the imbal- 
anced Fermi gas, where it leads to unphysical results for 
small chemical potential differences. In particularly, it 
predicts a negative polarization p = (n+ — n_ ) / (n+ — n_) 
for a positive chemical potential difference — fi~). 
However, by including Cooper-pair interaction effects 
with our RG, this problem is resolved. As a result, we 
also obtain very different results for the equation of state 
of an imbalanced Fermi gas compared to NSR theory 
and mean-field theory. The obtained pressure at nonzero 
temperatures can be measured in upcoming experiments. 
Moreover, the NSR theory is not accurate for determin- 
ing the ratio between the critical temperature and the 
chemical potential, for which it gives the same result as 
mean-field theory. Calculating self-energy effects with 
our RG theory for the balanced Fermi gas results in a 
critical temperature reduction of more than a factor of 
1.5, bringer it closer to Monte-Carlo results. Our method 
is complementary to a Fermi-liquid theory of the normal 
state, which is based on single-particle correlations. The 
Cooper-liquid theory calculates pair-correlation effects, 
and our method has the advantage that it can in principle 
be directly generalized to the superfluid state by perform- 
ing a RG for the pair condensed phase Moreover, it 
can also be extended to the mass-imbalanced case, which 
has a very rich phase diagram in the unitarity limit fo"l| . 
These are also projects for the future. 
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Appendix A: Useful relations 

Using contour integration, the Green's function of the 
Cooper pairs from Eq. Q can be written in the spectral 
form 

G A («w„,fc, fi A ) = - / dio ; , (Al) 

n J uj - icu n 

where lo + = u> + irj with rj 4. 0. The imaginary par t of the 
Green's function can be obtained analytically [52]. With 
the spectral representation, we can rewrite Matsubara 
sums over the pair Green's function as frequency inte- 
grals that are convenient for numerical evaluation. For 
example, we have 

^-^G A («w„,fc,^A) = 

n 

~ [ duiN B (u))lm[GA(u + ,k,tiA)}, (A2) 

where n is even and N-q(lj) = l/(e^ u — 1) is the bosonic 
distribution function. Moreover, the pair bubble diagram 
from Eq. ([9]) becomes 

U(k, n A ) = — ^2 Ga(«w„, k, ha) 2 = (A3) 

n 

^ J du}N B (u})Im[GA(^ + ,k, H A )}Re[G A {^ + ,k, (ia)}, 



where we used Eq. (|A1[) and the Kramers-Kronig relation 
to relate the real and imaginary part of the Cooper-pair 
Green's function. For the thermodynamic potential den- 
sity from Eq. (JT2J, we use that [Hf 



n 

i y^iVB(w)Im(log[-G A 1 ( W + ,fc, A/ A)]). (A4) 



Acknowledgements. — We thank Munekazu Horikoshi 
and Sylvain Nascimbene for kindly providing us with 
their experimental data, as well as Evgeni Burovski 
for sharing their Monte-Carlo data. This work is sup- 
ported by the Stichting voor Fundamcntcel Onderzoek 
der Materie (FOM) and the Nederlandse Organisatie 
voor Wetenschaplijk Onderzoek (NWO), as well as by 
the European Community's Seventh Framework Program 
ERC-2009-AdG under grant agreement 247142-MolChip. 



[1] M. Inguscio, W. Ketterle, and C. Salomon, eds., Ultra- 
cold Fermi Gases (IOS Press, Amsterdam, 2008), pro- 
ceedings of the International School of Physics Enrico 
Fermi, Course CLXIV, Varenna, 20-30 June 2006. 

[2] H. T. C. Stoof, K. B. Gubbels, and D. B. M. Dickerscheid, 
Ultracold Quantum Fields (Springer, Dordrecht, 2009). 

[3] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

[4] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 
Mod. Phys. 82, 1225 (2010). 

[5] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 
92, 040403 (2004). 

[6] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. 
Raupach, A. J. Kerman, and W. Ketterle, Phys. Rev. 
Lett. 92, 120403 (2004). 

[7] J. Kinast, S. L. Hemmer, M. E. Gehm, A. Turlapov, and 
J. E. Thomas, Phys. Rev. Lett. 92, 150402 (2004). 

[8] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, 
C. Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. 
Lett. 92, 203201 (2004). 

[9] T. Bourdel, L. Khaykovich, J. Cubizolles, J. Zhang, 
F. Chevy, M. Teichmann, L. Tarruell, S. J. J. M. F. 
Kokkelmans, and C. Salomon, Phys. Rev. Lett. 93, 
050401 (2004). 

[10] G. B. Partridge, K. E. Strecker, R. I. Kamar, M. W. Jack, 

and R. G. Hulet, Phys. Rev. Lett. 95, 020404 (2005). 
[11] T. L. Ho, Phys. Rev. Lett. 92, 090402 (2004). 
[12] S. Tan, Ann. Phys. 323, 2952 (2008). 



[13] L. Luo, B. Clancy, J. Joseph, J. Kinast, and J. E. 

Thomas, Phys. Rev. Lett. 98, 080402 (2007). 
[14] M. Horikoshi, S. Nakajima, M. Ueda, and T. Mukaiyama, 

Science 327, 442 (2010). 
[15] S. Nascimbene, N. Navon, K. J. Jiang, F. Chevy, and 

C. Salomon, Nature 463, 1057 (2010). 
[16] A. J. Leggett, Modern Trends in the Theory of Condensed 

Matter (Springer- Verlag, Berlin, 1980), p. 13. 
[17] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 

195 (1985). 

[18] C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, 

Phys. Rev. Lett. 71, 3202 (1993). 
[19] E. Burovski, N. Prokof'ev, B. Svistunov, and M. Troyer, 

Phys. Rev. Lett. 96, 160402 (2006). 
[20] H. Hu, P. D. Drummond, and X.-J. Liu, Nat. Phys. 3, 

469 (2007). 

[21] H. Hu, X.-J. Liu, and P. D. Drummond, New J. Phys. 

12, 063038 (2010). 
[22] M. W. Zwierlein, A. Schirotzek, C. H. Schunck, and 

W. Ketterle, Science 311, 492 (2006). 
[23] G. B. Partridge, W. Li, R. I. Kamar, Y. Liao, and R. G. 

Hulet, Science 311, 503 (2006). 
[24] Y. I. Shin, C. H. Schunck, A. Schirotzek, and W. Ketterle, 

Nature 451, 689 (2008). 
[25] K. B. Gubbels and H. T. C. Stoof, Phys. Rev. Lett. 100, 

140407 (2008). 

[26] K. B. Gubbels, M. W. J. Romans, and H. T. C. Stoof, 
Phys. Rev. Lett. 97, 210402 (2006). 



10 



[27] M. M. Parish, F. M. Marchetti, A. Lamacraft, and B. D. 
Simons, Nat. Phys. 3, 124 (2007). 
R. Stratonovich, Sov. Phys. Dok. 2, 416 (1958). 
J. Hubbard, Phys. Rev. Lett 3, 77 (1959). 

D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. 
Rev. Lett. 93, 090404 (2004). 

P. Pieri and G. C. Strinati, Phys. Rev. B 61, 15370 
(2000). 

P. Pieri and G. C. Strinati, Phys. Rev. B 71, 094520 
(2005). 

M. W. J. Romans and H. T. C. Stoof, Phys. Rev. A 74, 
053618 (2006). 

R. B. Diener, R. Sensarma, and M. Randeria, Phys. Rev. 
A 77, 023626 (2008). 

R. Haussmann, W. Rantner, S. Cerrito, and W. Zwerger, 
Phys. Rev. A 75, 023610 (2007). 
R. Haussmann, Phys. Rev. B 49, 12975 (1994). 

E. Burovski, E. Kozik, N. Prokof'ev, B. Svistunov, and 
M. Troyer, Phys. Rev. Lett. 101, 090402 (2008). 
J. Carlson, S. Y. Chang, V. R. Pandharipande, and K. E. 
Schmidt, Phys. Rev. Lett. 91, 050401 (2003). 
G. E. Astrakharchik, J. Boronat, J. Casulleras, and 
S. Giorgini, Phys. Rev. Lett. 93, 200404 (2004). 
A. Bulgac, J. E. Drat, and P. Magierski, Phys. Rev. A 
78, 023625 (2008). 



[41] S. Nascimbene, N. Navon, S. Pilati, F. Chevy, S. Giorgini, 
A. Georges, and C. Salomon, arXiv: 1012.4664 (2010). 

[42] K. G. Wilson and J. Kogut, Phys. Rep. 12, 75 (1974). 

[43] D. S. Fisher and P. C. Hohenberg, Phys. Rev. B 37, 4936 
(1988). 

[44] M. Bijlsma and H. T. C. Stoof, Phys. Rev. A 54, 5085 
(1996). 

[45] N. Hasselmann, S. Ledowski, and P. Kopietz, Phys. Rev. 

A 70, 063621 (2004). 
[46] J. P. Blaizot, R. Mendez-Galain, and N. Wschebor, Phys. 

Rev. E 74, 051117 (2006). 
[47] S. Floerchinger and C. Wetterich, Phys. Rev. A 77, 

053603 (2008). 

[48] Y. I. Shin, A. Schirotzek, C. H. Schunck, and W. Ketterle, 

Phys. Rev. Lett. 101, 070404 (2008). 
[49] M. C. Birse, B. Krippa, J. A. McGovern, and N. R. 

Walet, Phys. Lett. B 605, 287 (2005). 
[50] S. Diehl, H. Gies, J. M. Pawlowski, and C. Wetterich, 

Phys. Rev. A 76, 021602(R) (2007). 
[51] K. B. Gubbels, J. E. Baarsma, and H. T. C. Stoof, Phys. 

Rev. Lett. 103, 195301 (2009). 
[52] G. M. Falco and H. T. C. Stoof, Phys. Rev. A 75, 023612 

(2007). 



